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1.  INTRODUCTION 


A  data  structure  with  local  access  pr emits  an  access  pointer  Into  the 
structure  to  change  Its  value  from  one  location  only  to  an  adjacent  location 
in  the  structure.  For  example,  trees  are  typically  implemented  with  link 
fields  in  each  record  to  the  parent  and  the  children;  after  information  in 
a  record  is  retrieved,  the  next  record  accessed  is  one  of  these  adjacent 
records.  Even  arrays  may  have  local  access  if  implemented  by  orthogonal 
lists  [Kl]  rather  than  random  access  to  entries  by  address. 

Multihead  Turing  machines  model  storage  and  retrieval  operations  in  data 
structures  with  local  access:  tree  machines  for  trees,  multidimensional 
machines  for  arrays.  The  access  heads  of  the  machines  correspond  to  the  access 
pointers  into  the  data  structures .  At  each  step,  the  machine  can  shift  each 
access  head  only  to  an  adjacent  cell  in  its  worktape  storage.  Though  possibly 
unrealistic,  the  simple  input/output  and  control  conventions  are  not  essential 
for  our  results,  which  depend  only  on  the  storage  structure  and  the  local 
access  method. 

Paul,  Seiferas,  and  Simon  [PSS]  established  nonlinear  lower  bounds  on 
the  time  required  to  simulate  tree  machines  on-line  by  tree  machines  with 
fewer  access  heads  and  to  simulate  multidimensional  machines  on-line  by 
multidimensional  machines  with  fewer  access  heads.  In  essence,  increasing 
the  number  of  access  pointers  into  a  tree  or  an  array  can  save  time. 

We  present  efficient  simulations  of  multihead  machines  by  machines  with 
fewer  access  heads. 

Theorem  1.  Every  multihead  tree  machine  of  time  con^ laxity  t(n)  can  be 
simulated  on-line  by  a  tree  machine  with  two  access  heads  on  separate  work- 
tapes  in  time  O(t(n)log  t(n)/log  log  t(n>). 
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Theorem  2.  For  all  d  2  2  and  e  i  2,  every  multihead  e-dimensional  machine 
of  time  complexity  t(n)  can  be  simulated  on-line  by  a  d-dimeneional  machine 
with  two  access  heads  on  separate  worktapes  in  time  0(t(n)^  *  -  l/de^o^  t(n>). 

(Slightly  stronger  versions  of  Theorems  1  and  2  appear  as  Theorems  3  and 
4  below.) 

The  simulation  for  tree  machines  is  optimal  (within  a  constant  factor). 
Theorem  1  achieves  the  lower  bound  0 (t(n) log  t(n)/log  log  t(n))  of  Paul, 

Seif eras,  and  Simon. 

The  simulation  for  multidimensional  machines  is  not  known  to  be  optimal. 

The  known  lower  bound  on  the  time  required  for  simulating  e-dimensional 
machines  by  d-dimensional  machines  with  just  two  access  heads  on-line  is 
0(t(n)*  +  "  1^de  "  ®)  for  every  6  >  (d-l)/de(e+l)  and  d  i  2  ( PSS  ] .  Neverthe¬ 

less,  our  simulation  is  faster  than  the  simulation  of  Paul,  Seiferas,  and  Simon, 
which  for  e  *  d  at  2  runs  in  time  0(t(n)^ +  ”  s) ,  «  »  l/(d^(d-l) +d) . 

Furthermore,  our  simulation  employs  two  data  representations,  converting  from 
one  to  the  other  when  necessary.  This  technique  may  be  useful  in  other 
situations.  For  the  case  d  ■  1,  Stoss  [S]  generalized  a  classic  result  of 
Hennle  and  Stearns  [HS]  and  discovered  a  singulation  of  an  e-dimensional  machine 
by  a  one -dimensional  machine  with  two  access  heads  on  separate  worktapes  in 
time  0(t(n)2  ”  ^log  t(n)). 


] 
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2.  TREE  MACHINES 

Tree  machines  generalize  conventional  multihead  Turing  machines  with 
linear  tapes.  A  tree  worktape  is  a  collection  of  cells  organized  into  a 
complete  infinite  rooted  binary  tree.  The  root  has  two  children,  and  all 
other  cells  each  have  two  children  and  a  parent.  A  tree  machine  has  a 
finite-state  control,  a  read-only  linear  Input  tape,  a  write-only  linear 
output  tape,  and  a  finite  number  of  tree  worktapes,  each  of  which  has  a  finite 
number  of  access  heads.  At  each  step  the  machine  reads  the  symbols  in  the 
cells  on  which  the  input  and  worktape  heads  are  positioned,  writes  symbols 
from  a  finite  alphabet  on  these  worktape  cells  and  possibly  on  the  output 
tape  too,  and  possibly  shifts  each  head  to  an  adjacent  cell.  Initially,  all 
worktape  calls  hold  blanks,  and  every  worktape  head  is  positioned  on  the  root 
of  its  tape.  The  machine  can  detect  when  one  of  its  heads  is  at  the  root  and 
when  two  heads  scan  the  same  cell. 

Let  W  be  a  tree  worktape.  To  each  cell  of  W  assign  in  a  natural  way  a 
string  in  {0,1}*  called  the  location  of  the  cell.  Write  W(b)  for  the  cell 
of  W  at  location  b.  Let  X  denote  the  empty  string.  Then  W(X)  is  the  root 
of  W.  Write  W[b,r]  for  the  complete  subtree  of  W  of  height  r  rooted  at  W(b). 
Cell  W(a)  is  at  distance  r  from  cell  W(b)  if  an  access  head  requires  exactly 
r  steps  to  move  fromW(a)  to  W(b).  For  example,  the  leaves  of  W[b,r]  are 
at  distance  r  fromW(b). 

Let  T  be  a  tree  machine.  To  establish  Theorem  1,  we  devise  an  on-line 
simulation  of  T  by  a  tree  machine  S  with  two  heads  on  separate  worktapes. 
Recall  that  in  general,  machine  M'  simulates  machine  M  of  time  complexity 


4 


t(n)  on-line  In  tin*  f(t(n))  If  on  every  input  word,  if  *2  ...  S 

are  the  steps  at  which  M  reads  or  writes  a  symbol,  then  there  are 
corresponding  steps  ^  ^  •••  -  at  wbich  M>  read*  or  writes  the 

same  symbols,  and  every  s^  S  f(s^). 

Consider  a  confutation  of  N  steps  on  T.  For  simplicity,  we  present  a 
simulation  in  which  N  is  available  off-line.  The  simulation  can  be  converted 
to  an  on-line  simulation  with  time  loss  of  only  a  constant  factor  [PSS]. 

i 

(For  N*  -  1,2, 4, 8,...,  machine  S  repeats  the  simulation  from  the  beginning 
using  N'  for  the  length  of  the  computation  until  N'  is  at  least  as  large 
as  N.)  Without  loss  of  generality,  assume  that  every  worktape  head  of  T 
remains  within  distance  O(log  N)  of  the  root  of  its  tape.  Section  2.3 
justifies  this  assumption.  Also,  assume  that  the  h  access  heads  of  T  are  on 
separate  worktapes  W^, . . . .W^.  Section  2.2  handles  worktapes  with  multiple 
heads. 

2.1  Simulation 

Call  one  worktape  of  S  the  primary  tape  P,  the  other  the  secondary  tape. 
Section  2.2  replaces  the  secondary  tape  by  a  linear  tape.  Divide  the  cells 

of  P  Into  4h  tracks,  four  for  each  worktape  of  T,  numbered  . 

(worktape  contents),  . (head  positions),  C^,...,C^  (local  contents), 

(local  head  positions).  In  the  initial  configuration  of  S,  the 
blank  symbol  in  every  cell  is  Interpreted  as  a  blank  symbol  on  every  track. 

Set  r  -  (log  log  N)/2.  The  N  steps  of  the  computation  of  T  are  simulated 
in  N/r  phases  of  r  steps  each.  At  the  end  of  the  simulation  of  each  phase, 
track  Ct  of  P  holds  the  simulated  contents  of  W^  for  every  b,  track  Ct 
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of  P(b)  holds  Che  same  symbol  as  W^(b)  would.  Furthermore,  at  the  end  of 
each  simulated  phase,  track  indicates  a  path  from  the  root  of  P  to  P(x^) 
if  the  simulated  head  on  is  positioned  on 

To  simulate  one  phase,  S  first  determines  the  location  a^  of  the  ancestor 
at  distance  r  from  the  current  location  of  the  simulated  access  head  on  W^, 
setting,  a^  -  X  if  Wj^xp  is  within  distance  r  of  the  root  W^X).  Next,  S 
copies  the  contents  of  tracks  and  of  P[a^,2r]  to  the  corresponding  tracks 
and  of  P[X,2r],  using  the  secondary  tape  for  intermediate  storage. 

Observe  that  the  access  head  on  remains  in  U^[a^,2r]  during  this  phase  of 

r  steps.  Thus,  when  S  completes  these  copying  operations,  tracks  . 

of  P[X,2r]  hold  the  contents  of  all  cells  that  the  sinulated  heads  of  T 
could  access  during  the  phase. 

To  simulate  one  step  of  T,  machine  S  employs  track s  and  of  P[X,2r] 

to  determine  for  each  i  the  contents  of  the  cell  read  by  the  simulated  head 
on  U^.  Machine  S  then  records  the  new  contents  of  the  cell  on  track  and 
the  new  head  position  on  track  H^.  (The  management  of  track  is  entirely 
routine.)  At  the  end  of  the  phase,  for  each  i,  S  uses  the  secondary  tape 
to  copy  the  contents  of  tracks  and  of  P[X,2r]  back  onto  tracks 
and  H£  of  P[a.,2rl. 

Let  us  reckon  the  time  consumed  by  S .  At  the  beginning  and  end  of  each 
phase,  S  moves  its  head  on  the  primary  tape  to  each  P(a^)  (time  O(log  N) 
by  hypothesis  on  T).  It  copies  the  contents  of  h  trees  of  0(4r)  cells  each 
to  the  top  of  P  and  back  (time  0(4r)).  To  s inula te  one  step,  S  twice  moves 
the  head  to  a  cell  on  P[X,2r]  for  each  simulated  head  of  T  (time  0(r)). 
Therefore,  S  runs  in  time 

(N/r)[0(log  N)+0(4r)  +  r0(r)l  -  0(N  log  N/log  log  N). 
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2.2  Multihead  Tapes  and  a  Linear  Secondary  Worktape 

We  modify  S  to  simulate  T  when  T  has  multihead  tapes.  Clearly,  we  may 
assume  that  the  h  access  heads  of  T  are  on  the  same  worktape  W.  Divide  the 
primary  worktape  of  S  into  3h  +  l  tracks  called  C,  H^, . . .  ,H^,  C^,...,C^, 

holds  the  contents  of  the  simulated  worktape  W,  and  tracks 

specify  the  access  head  positions.  Tracks  C^,...,C^,  H| . are  used  to 

simulate  individual  steps. 

Machine  S  sets  r  *  (log  log  N)/2  and  simulates  T  in  N/r  phases  of  r 
steps  each.  To  simulate  a  phase,  S  determines  the  location  a^  of  the 
ancestor  at  distance  r  from  the  current  location  of  the  i-th  simulated  access 
head  of  T.  In  time  0(4r),  S  decides  for  each  pair  (i,j)  whether  W[a^,2r]  meets 
W[aj,2r].  If  so,  then  because  head  j  might  read  what  head  i  writes  during 
this  phase,  S  copies  track  C  of  both  Pta^r]  and  P[aj,2r]  to  the  same  C' 
track  of  P[X,2r] .  Since  the  time  to  simulate  one  step  once  the  copying 
operations  are  completed  is  0(hr)  *»  0(r),  this  modified  S  still  runs  in 
time  0(N  log  N/log  log  N) . 

(To  determine  whether  W[a^,2r]  meets  W[aj,2r],  S  marks  cells  of  P[aif2r] 
and  cells  of  P[aj,2r]  and  checks  whether  P[a^,2r]  meets  P[a^,2r],  In  its 
finite  state  control  S  constructs  a  graph  with  h  vertices  £l,...,h}  and  all 
edges  (i,j)  such  that  W[a^,2r]  meets  W[aj,2r].  The  connected  components  of 
this  graph  specify  which  P[a^,2r]  should  be  copied  to  the  same  C'  track  of 
P[X,2r].) 


,H^.  As  before,  at  the  end  of  the  simulation  of  a  phase,  track  C 
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Suppose  T  has  head-to-head  Jumps:  at  each  step  head  1  can  Jump  to  the 
cell  on  which  head  j  is  located.  In  terns  of  data  structures,  pointer  i 
can  be  assigned  the  current  value  of  pointer  j.  To  simulate  this  head-to- 
head  jump,  S  records  the  new  position  of  head  i  on  track  of  P[X,2hr] 
in  time  0(r).  When  completing  the  simulation  of  the  phase,  S  adjusts  the 
contents  of  track  in  time  O(log  N). 

We  can  replace  the  secondary  tree  worktape  of  S  by  a  linear 

(one-dimensional)  worktape  without  time  loss.  Initially,  S  writes  a  pattern 

of  length  0(2^7)  on  the  linear  tape  that  describes  a  depth-first  traversal  of 

a  complete  binary  tree  of  height  2r.  This  can  be  accomplished  by  using  the 

2r 

primary  tape  as  a  linear  tape  to  count  to  2  in  binary.  To  copy  the  contents 
of  a  tree  P[b,2r]  from  the  primary  worktape,  S  copies  the  contents  of  the  tree 
onto  this  pattern,  which  specifies  the  movement  of  the  head  on  the  primary 
worktape . 

We  summarize  these  alterations  to  S  in  a  form  that  subsumes  Theorem  1. 

Theorem  3.  Every  multihead  tree  machine  with  head-to-head  jumps  of 
time  complexity  t(n)  can  be  simulated  on-line  by  a  machine  with  a  tree 
worktape  and  a  linear  worktape,  each  with  one  access  head,  in  time 
O(t(n)log  t(n)/log  log  t(n)). 

2.3  Depth  Reduction 

We  establish  a  depth  reduction  lemma  for  multihead  tree  machines.  Paul 
and  Reischuk  [PR]  obtained  a  similar  result. 

Lemma  1.  Every  multihead  tree  machine  T  of  time  complexity  t(n)  can 
be  simulated  on-line  by  a  multihead  tree  machine  S  of  time  complexity  0(t(n)) 
whose  worktape  heads  r  main  with4  distance  0 (log  t(n))  of  the  root. 
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Proof.  As  usual,  assume  that  t  Is  a  power  of  2  and  is  available  off-line 
Also,  without  loss  of  generality,  assume  that  all  h  access  heads  of  T  are  on 
one  worktape  W.  The  simulator  S  has  several  access  heads.  Including  three 
for  each  head  of  T,  on  one  worktape  V.  Throughout  the  simulation  these 
heads  remain  within  distance  O(log  t)  of  the  root  of  V. 

Divide  V  into  two  tracks ,  numbered  0  and  1 .  In  the  initial  configuration 
a  blank  in  a  cell  of  V  is  interpreted  as  a  blank  on  both  tracks.  During 
the  simulation,  track  0  holds  the  contents  of  cells  of  W,  and  track  1  holds 
control  information. 

Set  L  ■  log(ht).  For  all  b,  call  W[b,L]  the  top  half  of  W[b,2L  +  l]  and 

the  rest  of  W[b,2L+l]  its  bottom  half.  Let  vo*vi»  *  *  *  *vht  136  t*le  l°cat*-on8 

of  cells  at  distance  L  from  the  root,  and  set  ■  V[vit2L  +  l]. 

Track  0  of  VQ  always  represents  the  cells  of  W[X,2L  + 1]  literally: 
for  every  cell  W(x)  in  W[X,2L  +  11,  the  symbols  in  W(x)  and  track  0  of 
V(VqX)  are  the  same.  (vqX  is  the  concatenation  of  the  binary  strings 
Vq  and  x.) 

For  j  >  0,  however,  only  the  bottom  half  of  holds  symbols  written 
in  W.  When  a  i*  X,  the  contents  of  W[a,2L  +  1]  correspond  to  V^.Vj,  where 
1  <  j,  if 

(i)  there  is  a  cell  V(Uj)  in  at  distance  L  from  V(v^)  such  that 
Vj  is  written  on  track  1  of  V[u^,L]; 

(ii)  Uj  is  written  on  track  1  of  V[Vj,L];  and 

(iii)  for  all  x,  if  W(ax)  is  in  the  top  half  of  W[a,2L  +  l],  then  the 
symbol  in  W(ax)  is  the  same  as  the  symbol  on  track  0  of  V(u^x), 
whereas  if  W(ax)  is  in  the  bottom  half  of  W[a,2L+i],  then  the 
symbol  in  W(ax)  is  the  same  as  the  symbol  on  track  0  of  V(v^x). 
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Throughout  the  simulation,  S  maintains  the  contents  of  W  through  these 
correspondences.  Suppose  either  that  W[a,2L+l]  corresponds  to  or 

that  a  »  X  and  j  ■  0.  If  W(b)  ■  W(ac)  is  at  distance  L  +  1  from  W(a)  in 
W[a,2L  +  l]  and  the  bottom  half  of  W[b,2L  +  l]  is  nonblank,  then  track  1  of 
V[VjC,L]  holds  some  v^,  and  W[b,2L  +  l]  corresponds  to  See  Figure  1. 

Observe  that  the  relative  position  of  v^c  in  Vj  is  the  same  as  the  relative 
position  of  b  in  W[a,2L  +  l].  Conversely,  if  V(VjC)  is  at  distance  L+l 
from  V(Vj)  in  and  track  1  of  V[VjC,L]  holds  some  v^,  then  W[ac,2L+l] 
corresponds  to  ,V^. 

For  every  access  head  H  of  T,  machine  S  has  access  heads  A  (the  top 

m  m 

head)  and  Bm  (the  bottom  head).  Suppose  tf[a,2L+l]  corresponds  to 

To  simulate  Hm  reading  cell  W(ax)  using  V^,Vj,  the  top  head  Affl  is  at  cell 

V(UjX)  and  the  bottom  head  at  V(VjX).  If  W(ax)  is  in  the  top  half  of 

W[a,2L  +  l],  then  the  top  head  Am  is  in  V[u^  ,L] ,  and  S  retrieves  the  symbol 

in  W(ax)  from  track  0  of  V(UjX).  If  W(ax)  is  in  the  bottom  half  of  W[a,2L  +  l], 

then  S  uses  the  track  0  contents  of  V(v^x)  read  by  the  bottom  head  Bffl.  A 

third  access  head  Cffl  is  used  by  S  as  a  unary  counter  to  determine  whether 

W(ax)  is  in  the  top  half  or  the  bottom  half  of  W[a,2L+l]. 

In  the  step-by-step  simulation  of  T,  there  are  two  situations  in 

which  A  and  B  require  reorientation, 
m  m  '  - 

(1)  Let  W[a,2L  +  l]  correspond  to  V,,V.,  and  let  A  ,B  use  V,,V.. 

i  j  m  m  i  j 

Suppose  attempts  to  visit  a  child  of  a  leaf  W(az)  of  W[a,2L  +  l].  Let 
W (b)  be  the  ancestor  of  W(a z)  at  distance  L  from  W(az).  Similarly,  let 


V(VjC)  be  the  ancestor  of  V(VjZ)  at  distance  L  from  V(VjZ).  Define  d  so 
that  cd  ■  z.  See  Figure  1.  If  Crack  1  of  V[VjC,L]  contains  some  v^. 
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then  W{b,2L  +  l]  corresponds  to  Vj.V^.  In  this  case  S  sends  Am  from  V(UjZ) 

to  V(VjZ)  In  Vj,  sends  Bffl  from  V(VjZ)  to  V(vkd)  in  V^,  and  resumes 

the  simulation  with  A  ,B  using  V,,V.  .  Otherwise,  if  track  1  of  V(v,c,Ll 

mm  J  *  j 

is  blank,  then  S  marks  the  first  unmarked  V(vk),  writes  vfe  in  track  1  of 
V[VjC,L]  writes  v^c  in  track  1  of  V[vk,L],  and  proceeds  as  above,  for  now 
W[b,2L+l]  corresponds  to  Vj,Vk.  (The  contents  of  track  0  of  are 
completely  blank.) 

(2)  Let  W(b,2L  +  ll  correspond  to  Vj,Vk,  and  let  Am,Bm  use  V^.V^. 
Suppose  that  Hm  attempts  to  visit  the  parent  of  W(b).  Let  W(a)  be  the 
ancestor  of  W(b)  at  distance  L+l  from  W(b).  Define  c  so  that  b  -  ac,  and 
set  u^  *  VjC.  From  track  1  of  V[Vj,L],  S  retrieves  u^,  where  V(Uj)  is  in 
some  .  We  know  that  W(a,2L+l]  corresponds  to  V^,V^.  S  sends  head  B^ 
from  V(vfc)  to  VOi^)  and  head  Affl  from  VO^)  to  V(UjC),  and  resumes  the 
simulation  with  A  ,B„  using  V. ,V4. 

mm  i  j 

Since  for  each  access  head  H  of  I,  machine  S  simulates  at  least 

m 

L  *  ®(log  t)  steps  between  reorientations  of  A^.B^,  and  each  reorientation 
takes  time  0(log  t),  we  conclude  that  S  simulates  T  in  linear  time.  C 

Corollary.  Every  multihead  tree  machine  with  head-to-head  jumps  of 
time  complexity  t(n)  can  be  simulated  on-line  by  a  tree  machine  with 
head-to-head  jumps  of  time  complexity  0(t(n))  whose  worktape  heads  remain 
within  distance  0(log  t(n))  of  the  root. 


Proof.  We  modify  the  simulator  S  described  In  the  proof  of  Lemma  1. 

Assume  that  the  tree  machine  T  has  h  heads  on  one  worktape  W.  Both  S  and  T 

permit  head-to-head  jumps.  Set  L  ■  log(ht). 

For  each  worktape  head  H  of  T,  machine  S  has  several  heads  on  the  same 

n 

worktape  V,  Including  A'  (ancestor  top  head)  and  B1  (descendant  bottom  head) 

m  m 

In  addition  to  the  top  head  A  and  the  bottom  head  B_  described  earlier. 

Ill  in 

The  heads  A*  and  B'  will  assist  In  reorientations.  We  now  specify  when  A' 
o  n  m 

and  B'  are  ready, 
m  — — i 

Suppose  A  ,  B  use  V.,  V.  to  simulate  H  In  W[b,2L+l],  cell  W(a) 
mm  j  k  m 

Is  the  ancestor  of  W(b)  at  distance  L+l  from  W(b),  and  b  *  ac:  the  top  head 

A  is  at  V(v.cx),  and  the  bottom  head  B  Is  at  V(v,x).  See  Figure  1.  If  A' 
m  J  m  k  m 

is  ready,  then  A'  is  at  V(u.c),  a  descendant  of  V(v,),  such  that  u,  is 
m  j  1  J 

written  on  track  1  of  V[Vj,L].  Note  that  If  A^,  Bm  used  V^,  Vj  to  simulate 
Hm  at  W(b)  in  W[a,2L  +  1],  then  Am  would  be  at  V(UjC). 

Now  suppose  further  that  W(be)  Is  at  distance  L+l  from  W(b)  and  Bm  is 
at  V(vfeex)  in  the  bottom  half  of  Vfc.  If  B^  is  ready,  then  Is  at  V(v^x) 

In  V^,  where  W[be,2L+l]  corresponds  to  V^,  V^.  Note  that  if  Am,  Bm  used 

Vk.  V,  to  simulate  Hm  in  W(be,2L+l],  then  Bm  would  be  at  V(v^x). 

The  on-line  simulation  of  T  by  S  proceeds  as  before.  Let  W[b,2L+l] 
correspond  to  V.,  V,  ,  and  let  A  ,  B  use  V.,  V.  .  When  H  shifts  to  an 

J  K  min  J  K  m 

adjacent  cell  in  W[b,2L  +  l],  heads  A  and  B  shift  to  the  corresponding 

cells  in  V.  and  V.  .  Furthermore,  if  B'  in  V.  is  ready,  and  B  remains  in 
j  k  m  i  m 

the  bottom  half  of  V.  ,  then  B'  also  shifts  to  an  adjacent  cell  in  V.  and 

lc  m  * 


remains  ready. 
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To  simulate  a  jump  of  H  to  H  ,  all  heads  of  S  for  H  jump  to  the 

IB  p  A 

corresponding  heads  for  H  :  head  A  jumps  to  A  ,  head  B_  to  B_ ,  head  A' 

P  q  p  01  p  A 

to  A',  and  head  B'  to  B'.  If  A*  is  ready,  then  S  marks  A*  ready;  if  B' 

p  A  p  p  A  p 

is  ready,  then  S  marks  B*  ready. 

A 

Suppose  H  shifts  to  the  parent  of  W(b).  If  A'  is  ready  at  a  descendant 

A  A 

of  v(v. ),  then  head  B  jumps  to  A.  head  A  jumps  to  A',  and  the  simulation 

X  A  A  A  A 

continues  with  A  ,  B  using  V.,  V.;  after  this  reorientation,  heads  A' 

A  A  X  j  A 

and  B^  are  no  longer  ready.  If  A^  is  not  ready  when  Hm  shifts  to  the  parent 

of  W(b),  then  a  widespread  reorientation  is  necessary  to  make  Am>  Bm  use 

Vt,  Vj.  Widespread  reorientations  are  discussed  below. 

Suppose  that  W(be)  is  at  distance  L+l  from  W(b)  and  H  shifts  to  a  child 

A 

of  a  leaf  W(bex)  of  W[b,2L+l].  If  B*  is  ready  in  V,,  then  head  A  jumps 

m  x  m 

to  B  ,  head  B  jumps  to  B1,  and  the  simulation  continues  with  A  ,  B  using 
m  m  m  mm 

Vk.  ;  after  this  reorientation,  heads  Aa  and  Ba  are  no  longer  ready. 

Suppose  Ba  is  not  ready  when  Hm  shifts  to  a  child  of  W(bex).  If  W[be,2L+l] 

corresponds  to  V^,  for  some  V^,  then  a  widespread  reorientation  makes 

A  ,  B  use  V,  ,  V,.  Otherwise,  the  bottom  half  of  W[be,2L+l]  has  not  been 
mm  k  x 

visited  previously,  and  an  initializing  reorientation  is  necessary  to 

initialize  a  new  V.  to  make  A  .  B_  use  Vv,  V.. 

x  mm  k  x 

As  in  the  proof  of  Leona  1,  for  an  initializing  reorientation  S  marks 
the  first  unmarked  V(v^),  writes  Vj  in  track  1  of  V[vke,L],  and  writes  vfce 
in  track  1  of  V[Vj,L].  Then  S  sends  head  A^  to  V(VjCe),  head  Aa  to  V(vkex) 
in  V,  ,  and  head  B  to  V(v.x)  in  V,;  it  marks  A'  ready  and  B'  not  ready. 
Furthermore,  for  all  other  H^,  it  makes  ancestor  top  head  A^  ready  and. 
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when  possible,  descendant  bottom  head  ready.  Observe  that  L  cells  in 
the  top  half  of  W[be,2L+l]  oust  have  been  visited  previously  because  the 
cells  visited  by  heads  of  T  always  form  a  connected  region  of  W.  Consequently, 
during  the  simulation  of  T  by  S,  there  are  only  0(t/L)  initializing 
reorientations,  each  of  which  takes  time  0(L).  The  total  time  for  initializing 
reorientations  is  0(t). 

When  H_  induces  a  widespread  reorientation,  S  uses  information  in  track  1 
m 

to  send  A  ,  B  to  the  appropriate  cells  in  subtrees  V . ,  V.  or  V,  ,  V.. 

mm  1  j  k  i 

Moreover,  for  every  p  *  l,...,h,  machine  S  makes  ancestor  top  head  A* 

P 

ready  and,  when  possible,  descendant  bottom  head  ready.  It  is  crucial 
to  realize  that  S  can  simulate  at  least  L  steps  of  T  between  widespread 
reorientations.  For  example,  suppose  that  inznediately  after  a  widespread 
reorientation,  A  ,  B  use  V.,  V.  to  simulate  H  lnW[b,2L+l],  and  B1  is 

p  p  j  K  p  p 

not  ready.  Then  either  Bp  is  in  the  top  half  of  V^,  or  Bp  is  in  V[vke,L] 
in  the  bottom  half  of  for  some  W(be)  at  distance  L  +  1  from  W(b),  and  there 
is  no  Vj  such  that  W[be,2L+l]  corresponds  to  V^,  V^.  In  the  former  case, 

Hp  is  at  distance  at  least  L+l  from  a  leaf  of  W[b,2L+l].  In  the  latter 
case,  if  Hp  visited  a  leaf  of  W[b,2L+l],  then  it  would  induce  an  initializing 
reorientation,  not  a  widespread  reorientation;  after  the  initializing 
reorientation,  all  ancestor  top  heads  and  as  many  descendant  bottom  heads  as 
possible  would  be  ready  again. 

Therefore,  since  S  simulates  at  least  L  steps  of  T  between  widespread 
reorientations,  there  are  only  0(t/L)  widespread  reorientations,  each  of 
which  takes  time  0(L).  The  total  time  for  widespread  reorelntations  is 
0(t) ,  and  S  simulates  T  in  linear  time.  □ 
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3.  MULTIDIMENSIONAL  MACHINES 

I 

The  Cwo  tape  simulation  of  a  multidimensional  machine  devised  by  Paul, 
Seiferas,  and  Simon  [PSS]  uses  pages  of  filed  size  and  a  hierarchy  of  simula¬ 
tion  procedures  to  maintain  the  locations  of  the  pages.  To  develop  a  faster 
simulation,  we  employ  a  directory  and  a  mechanism  for  handling  pages  of 

* 

variable  size  [L].  In  addition,  our  simulator  uses  both  trie  representations 
and  literal  representations  for  contents  of  the  multihead  machine's  vorfctapes. 
During  the  simulation,  it  may  convert  the  contents  of  a  page  from  one 
representation  into  the  other. 

To  each  cell  of  a  d -dimensional  worktape  assign  in  the  usual  way  a 
d-tuple  of  integers  called  the  coordinates  of  the  cell.  The  origin  is  the 
cell  whose  coordinates  are  all  zero.  A  d -dimensional  Turing  machine  has  a 
finite-state  control,  a  read-only  linear  input  tape,  a  write-only  linear 
output  tape,  and  a  finite  number  of  d-dlmensional  worktapes,  each  of  which 
has  a  finite  number  of  access  heads.  At  each  step  the  machine  reads  the 
symbols  in  the  cells  on  which  the  input  and  worktape  heads  are  positioned, 
writes  symbols  from  a  finite  alphabet  on  these  worktape  cells  and  possibly 
on  the  output  tape  too,  and  possibly  shifts  each  head  to  an  adjacent  cell. 
Initially,  all  worktape  cells  hold  blanks,  and  every  worktape  head  is 
positioned  on  the  origin  of  its  tape. 

On  a  d-dimensional  worktape,  a  box  is  a  set  of  cells  that  form  a 
d -dimensional  cube.  The  volume  of  a  box  is  the  number  of  cells  in  it. 

The  base  cell  of  a  box  is  the  cell  whose  coordinates  are  the  smallest. 

The  relative  position  of  a  cell  X  in  a  box  B  is  the  list  of  coordinates 


of  X  when  the  base  cell  of  B  is  taken  as  the  origin. 
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Fix  integer*  d  2  2,  e  *  2,  and  h.  Let  E  be  an  e -dimensional  machine 
with  h  heads  on  one  worktape.  To  establish  Theorem  2,  we  design  a 
d-dimsnsional  simulator  D  with  one  head  on  a  d-dimensional  worktape  and 
one  head  on  a  linear  worktape.  Technical  details  omitted  from  the 
presentation  of  the  simulation  in  Section  3.1  appear  in  Sections  3.2 
and  3.3. 

3.1  Simulation 

Consider  a  computation  of  N  steps  on  E.  As  before,  assume  that  N  is 
available  off-line.  The  coordinates  of  every  cell  accessed  by  E  during  the 
computation  can  be  written  as  a  binary  string  of  length  O(log  N) .  For 
convenience  write  log^N  for  (log  N)^. 

Define  the  function  tt  by 

n(x)  .  2ri°*  x1  ; 

if  x  is  not  a  power  of  2,  then  tt  maps  x  to  the  next  larger  power  of  2.  Set 

/ 

,1/de 


{ 


(N  log  N)‘ 


(M/log  N) 


1/4 


if  d  i  3  or  e  i  3 

if  d  -  e  -  2 


t  -  TT«b  log  N)1/d), 
u  -  TT(b*^d), 

u*  -  4  TT((3*h  N  log  N)l/d)  . 

At  the  cost  of  Introducing  another  constant  factor,  we  may  assume  that  b 
is  an  integer. 


v 


Partition  the  vorktape  of  &  into  pairwise  disjoint  boxes  of  side  b 
called  blocks .  The  position  of  a  block  is  the  list  of  coordinates  of  its 
base  cell. 

Divide  every  cell  of  the  d-dimensional  vorktape  of  D  into  three 
tracks  numbered  0,1,2.  In  the  initial  configuration  of  D,  the  blank  symbol 
in  every  cell  is  interpreted  as  a  blank  on  each  track. 

A  page  is  a  box  on  the  d-dlmensional  vorktape  whose  side  is  a  power  of 
2  between  t  and  u  that  has  a  nonblank  cell  counter  on  track  0.  We  describe 
how  the  contents  of  a  page  P  represent  the  contents  of  a  s inula ted  block  B. 
(Abusing  terminology,  we  shall  say  that  P  represents  B  when  the  configurations 
are  implicitly  specified.)  First,  the  value  of  the  nonblank  cell  counter 
must  be  at  least  the  number  of  nonblank  cells  in  B.  To  maintain  the  contents 
of  these  cells,  P  uses  either  a  binary  Internal  trie  on  track  1  or  a 
literal  representation  on  track  2. 

An  internal  trie  (sometimes  called  a  digital  trie  [K3,  p.  489])  is  a 
collection  of  nodes  organized  into  a  binary  tree  according  to  the  value 
of  a  binary  key:  the  location  of  each  node  v  in  the  tree,  written  as  a 
binary  string  0(v)  In  a  natural  way  as  in  Section  2,  is  an  initial  segment 
of  the  key.  In  the  page  P  a  node  is  a  box  of  fixed  volume  0(log  N)  together 
with  its  contents.  In  the  internal  trie  representation  of  B  there  is  a  node 
for  every  nonblank  cell  X  in  B,  and  the  coordinates  of  cells  are  used  as 
keys.  The  node  v  for  cell  X  holds  on  track  1 

the  coordinates  of  X  (the  key,  written  as  a  binary  string), 
the  contents  of  X,  and 

the  relative  positions  of  the  left  and  right  children  of  v  (if  any)  in  P. 
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The  number  of  nodes  equals  the  number  of  nonbleak  cells  in  B.  (We  may 
assume  that  E  never  writes  a  blank  on  a  vorktape  cell.)  See  Figure  2. 

The  internal  trie  representation  Is  particularly  useful  when  B  has  few 
nonblank  cells . 

The  literal  representation  Is  straightforward.  For  every  cell  X  of  B 
there  Is  a  cell  Y  of  P  whose  relative  position  In  P  is  easily  calculated 
from  the  coordinates  of  X  such  that  Y  holds  the  same  symbol  (on  track  2) 
as  X  does.  The  literal  representation  is  most  efficient  when  B  has  many 
nonblank  cells. 

On  the  d-dimensional  worktape  D  has  a  mass  store,  a  directory,  and  a 
working  area.  The  mass  store  is  a  box  of  side  u*  that  holds  the  pages. 

The  directory  is  a  box  of  side  u*/log  N  whose  contents  maintain  the  positions 
of  pages  in  the  mass  store.  During  the  simulation,  D  uses  the  position  of 
a  block  B  in  E  to  retrieve  from  the  directory  the  position  of  a  page  that 
represents  B.  (The  directory  is  described  in  Section  3.2.)  The  working 
area  comprises  3^  contiguous  boxes  of  side  u  at  fixed  positions;  call 
these  boxes  slots.  For  1  »  l,...^^,  the  position  of  the  i-th  slot  can 
be  calculated  quickly.  The  coordinates  of  every  cell  in  the  mass  store,  the 
directory,  and  the  working  area  can  be  specified  by  a  string  of  length 
O(log  N). 

The  linear  worktape  of  D  is  used  for  routine  arithmetic  calculations, 
for  counting,  for  maintaining  the  coordinates  of  the  simulated  heads  of  E, 
and  for  copying  pages  between  the  mass  store  and  the  working  area. 
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Machine  D  simulates  N  steps  o£  E  in  N/b  phases  of  b  steps  each.  At 
the  beginning  of  the  simulation  of  each  phase,  D  determines  the  positions 
of  the  at  most  3%  blocks  that  surround  the  h  simulated  worktape  heads  of  E ; 
during  the  phase  the  heads  remain  within  these  blocks.  For  each  of  these 
blocks  Bj ,  D  calls  procedure  LOCATEPAGE,  which  uses  the  directory  to  locate 
the  page  that  represents  Bj .  (LOCATEPAGE  is  described  in  Section  3.2.) 

Next,  D  copies  the  contents  of  P^  into  the  j-th  slot  in  the  working  area. 

In  its  finite-state  control  D  maintains  the  correspondence  between  blocks 
and  slots;  on  the  linear  worktape  it  records  the  positions  of  the  corresponding 
pages  Pj  in  the  mass  store.  At  the  end  of  the  phase,  D  copies  the  contents 
of  pages  from  the  working  area  to  the  mass  store:  D  copies  back  to  just 
the  box  in  that  was  written  at  the  beginning  of  the  phase.  To  send  the 
head  on  the  d-dimensional  worktape  to  the  proper  locations,  D  uses  the  linear 
worktape  as  a  counter. 

During  a  phase,  to  slnulate  one  step,  D  first  decides  for  each  simulated 
head  of  E  the  block  in  which  lies.  Next,  D  recalls  which  slot 
in  the  working  area  corresponds  to  B^.  Using  the  contents  of  Q^,  machine  D 
determines  the  contents  of  the  cell  X^  read  by  and  the  new  symbol  chat 

writes  on  X^.  If  represents  B^  literally,  then  D  calculates  the  relative 
position  of  the  cell  Y^  in  that  represents  X^,  reads  the  contents  of  Y^ 

(on  track  2),  and  writes  the  new  symbol  on  track  2  of  Y^.  If  employs  an 
internal  trie  to  represent  B^,  then  D  uses  the  coordinates  of  X^  to  find  the 
node  that  holds  the  contents  of  X^  on  track  1  and  records  the  new  contents 
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of  X^  la  this  node.  If  X^  has  not  bean  visited  previously,  then  It  holds  a 

blank,  and  D  creates  a  node  for  X^  at  an  appropriate  location  In  the  trie: 

D  adds  a  node  v  to  the  Internal  trie  for  which  P(v)  is  the  shortest  initial 

segment  of  the  coordinates  of  X^,  written  In  binary,  that  is  not  the  binary 

string  0(v')  for  a  node  v'  already  in  the  trie. 

Let  us  assess  the  time  that  D  takes  to  simulate  one  step.  First,  D 

2 

determines  the  coordinates  of  the  worktape  cells  read  by  E  (time  O(log  N) 

for  arithmetic  calculations  on  the  linear  worktape  because  this  worktape  has 

just  one  head).  If  page  Q  in  the  working  area  uses  a  literal  representation, 

then  to  retrieve  or  to  write  the  contents  of  a  cell,  D  computes  the  relative 

2 

position  of  the  representing  cell  in  Q  (time  O(log  N))  and  moves  its  head 

2 

across  the  working  area  (time  0(u))  for  a  total  of  0 (u  + log  N)  time.  If  Q 
uses  an  internal  trie  representation,  then  the  head  on  the  d-dimensional 
worktape  visits  O(log  N)  nodes  of  Q  -  at  most  one  node  for  each  bit  of  the 
coordinates  of  the  cell.  For  each  node,  D  copies  its  contents  onto  the  linear 
worktape  (time  O(log  N) )  and  moves  its  head  across  the  working  area  (time 
0(u))  for  a  total  of  O((u+log  N)log  N)  time.  Thus,  each  phase  takes  time 
0  (b  (u  +  log  N))  to  simulate  b  steps 

+  O(log^N)  to  calculate  the  positions  of  the  3%  blocks  that 
surround  the  heads  of  E 

+  0  (u*)  to  move  the  head  across  the  mass  store  to  each  of  3%  pages 
+  0(ud)  to  copy  the  contents  of  3eh  pages  of  volume  ud  between 
the  mass  store  and  the  working  area 
+  the  time  taken  by  calls  to  LOCATE PAGE . 

In  Section  3.2  we  show  that  during  the  simulation,  the  total  time  consumed 
by  calls  to  LOCATE PAGE  is  0((N/b)u*)  +  0(Nu  log  N).  Therefore,  the 
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simulation  runs  In  time 


(N/b)[0(u*)  +  0(u  )  +  0(bu  log  N)]  +  0(Nu  log  N) 

■  0(Nu*/b)  by  definition  of  b 

-  0(N1+  1/d  "  1/de  log1/d  “  l/de  N)  If  d  a  3  or  e  *  3 


KN5'4 


log3/4  N)  if  d  -  e  -  2, 


The  simulation  is  routinely  modified  to  handle  head-to-head  jumps  by  E. 


Theorem  4.  For  all  d  i  2  and  e  2  2,  every  oultihead  e-dimenslonal  machine 
with  head-to-head  jumps  of  time  complexity  t(n)  can  be  sinulated  on-line  by 
a  machine  with  a  d-dimensional  worktape  and  a  linear  worktape,  each  with  one 
access  head,  in  time 

0(t(n)1  +  1/d  ‘1/delog1/d  -1/det(n))  if  d  *  3  or  e  *  3, 
O(t(n)5/4log3/4t(n»  if  d  -  e  -  2. 

3.2  IX)  CATE  PAGE  and  the  Directory 

The  directory  of  D  is  a  box  of  side  u*/log  N  on  the  d -dimensional  worktape. 
To  maintain  the  positions  of  pages  whose  contents  represent  contents  of  blocks, 
the  contents  of  the  directory  are  organized  into  an  internal  trie  that  uses 
the  positions  of  blocks  as  keys.  For  every  nonblank  block  B  there  is  a 
node  v  in  the  directory  of  volume  0(log  N)  that  holds 
the  position  of  B, 

the  position  of  the  page  P  that  represents  B,  and 
the  positions  of  the  left  and  right  children  of  v  (if  any) 
such  that  0(v)  is  an  Initial  segment  of  the  position  of  B  (written  as  a  binary 
string  of  length  0(log  n)).  We  say  that  P  is  assigned  to  B.  A  page  in  the 
mass  store  is  active  if  it  is  assigned  to  some  block. 


l 


I 

£ 
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Throughout  the  simulation,  several  relationships  among  page  sizes  and 
nonblank  cell  counter  values  remain  invariant.  We  assert  these  relationships 
as  a  lesma. 

Leassa  2.  During  the  simulation  of  E  by  D,  the  following  always  hold. 

(a)  If  the  nonblank  cell  counter  of  a  page  P  has  value  m,  then  P 
has  side  min{u,TT((m  log  N)^d)). 

(b)  Let  page  P  represent  block  B.  If  P  has  side  u,  then  P  represents 
B  literally.  If  the  side  of  P  is  smaller  than  u,  then  P  uses  an  internal 
trie. 

(c)  There  are  at  most  3*h  N/b  active  pages. 

(d)  The  sum  of  the  nonblank  cell  counter  values  of  the  active  pages  is 
at  most  3*h  N. 

d  g 

Note  that  since  u  *  b  ,  if  page  P  has  side  u,  then  P  has  at  least  as 
many  cells  as  a  block  B  and  can  contain  a  representing  cell  for  every  cell  of 
B  in  a  literal  representation.  If  P  has  side  rr((m  log  N)^d) ,  where  m  is 
the  value  of  its  nonblank  cell  counter,  then  the  at  most  m  nodes  of  volume 
O(log  N)  in  the  internal  trie  representation,  one  for  each  nonblank  cell 
of  B,  fit  in  P. 

Furthermore,  by  Lemma  2(c),  the  directory  can  accomodate  the  nodes 
for  the  active  pages.  For  each  active  page  there  is  a  node  of  volume 
O(log  N).  The  implicit  constant  can  be  chosen  so  that 

(3*h  N/b)0(log  N)  £  (u*/log  N)d. 


22 


Now  we  describe  procedure  LOCATEPAGE.  Given  Che  position  of  a 
block  B,  LOCATEPAGE  employs  Che  directory  to  find  Che  page  P  in  Che  mass 
store  Chat  represents  B.  If  no  page  is  assigned  Co  B ,  then  LOCATEPAGE 
calls  procedure  ALLOCATE  to  produce  a  blank  box  of  side  t  in  Che  mass  store 
and  sets  Che  value  of  Che  nonblank  cell  counter  of  P  to  b;  since  B  must 
be  completely  blank  (it  has  not  been  previously  visited),  P  already 
represents  B. 

If  LOCATEPAGE  found  P  in  Che  mass  score  without  calling  ALLOCATE ,  then 
LOCATEPAGE  ensures  Chat  P  is  large  enough  to  include  all  nodes  that  D  may 
add  to  che  internal  crie  of  P  during  the  phase.  Let  p  be  the  side  of  P  and 
m  be  che  value  of  its  nonblank  cell  councer.  LOCATEPAGE  secs  che  value  of 
the  nonblank  cell  counter  of  P  to  m'  *  mintm+  b,be} .  Let 
p'  ■  mLnCutTT^m' log  N)^d)}.  If  p  <  p',  then  the  contents  of  P  are  copied 
into  a  new  blank  box  of  side  p'  (found  in  the  mass  store  by  ALLOCATE)  to 
produce  a  larger  page  P*  that  represetns  B  with  an  internal  trie.  If 
p  <  p'  -  u,  then  the  Internal  trie  representation  on  track  1  of  P'  is 
converted  into  a  literal  representation  on  track  2.  From  every  node  of  the 
internal  trie,  LOCATEPAGE  extracts  the  coordinates  of  a  cell  X  of  B  and 
copies  the  contents  of  X  from  track  1  into  track  2  of  the  cell  of  P1  that 
represents  X  in  the  literal  representation.  The  other  cells  of  B  are  blank, 
and  track  2  of  the  corresponding  cells  of  P'  remains  blank.  Finally, 
LOCATEPAGE  modifies  the  directory  to  assign  P'  to  B.  It  is  straightforward 
to  confirm  that  the  changes  made  by  LOCATEPAGE  satisfy  the  assertions  of 
Lemma  2. 


Let  us  evaluate  the  time  taken  by  calls  to  L0CAXEPAGE.  Section  3.3 

2 

demonstrates  that  ALLOCATE  operates  In  time  O(log  N).  In  the  directory, 
to  move  the  head  from  one  node  to  another  takes  time  proportional  to 
u*/log  N,  the  side  of  the  directory.  Consequently,  to  determine  the 
position  of  a  page  assigned  to  a  block  takes  time 

O(log  N)[0(log  N)  +  O(u*/log  N) ]  -  0(u*) 

because  O(log  N)  nodes,  each  of  volume  O(log  N)  are  accessed.  The  arithmetic 

2 

calculations  performed  by  LOCATEPAGE  rake  time  O(log  N). 

We  assess  the  time  for  converting  internal  trie  representations  into 
literal  representations  separately.  Let  page  P  represent  block  B,  and 
suppose  D  converts  the  internal  representation  in  P  into  a  literal  representa¬ 
tion.  Let  m  be  the  value  of  the  nonblank  cell  counter  of  P  when  D  performs 
the  conversion.  According  to  the  definition  of  LOCATEPAGE,  m  -  9(ud/log  N) . 
For  each  of  at  most  m  nonblank  cells  of  B,  machine  D  accesses  O(log  N)  nodes 
of  P  in  time  O((u+log  N)log  N)  .  By  Lemma  2(d),  D  converts  at  most 
0(N/m)  ■  0((N  log  N)/ud)  pages  during  the  simulation.  Therefore,  the  total 
time  for  converting  representations  is 

0((N  log  N)/ud)0(ud/log  N)0((u  +  log  N)log  N)  -  0(Nu  log  N) . 

We  conclude  that  the  total  time  consumed  by  calls  to  LOCATEPAGE  during 
the  simulation  is 

0(N/b)[0(u*) +0(log2N)l +0(Nu  log  N)  -  0(Nu*/b) +0(Nu  log  N). 
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3.3  Storage  Allocation 

Machine  D  keeps  a  free  storage  list,  a  list  of  addresses  of  blank 
boxes  in  the  ness  store.  For  q  -  1,2,4, ... ,u*/2,  u*,  the  free  storage  list 
has  positions  of  at  most  2**-l  boxes  of  side  q  in  the  mass  store.  Initially, 
the  free  storage  list  holds  the  position  of  the  mass  store,  a  single  box 
of  side  u*. 

Procedure  ALLOCATE  employs  a  buddy  system  [Kl]  to  allocate  a  blank  box 
with  a  desired  side  in  the  mass  store.  To  obtain  a  blank  box  of  side  r, 
a  '>ower  of  2,  ALLOCATE  finds  the  position  of  a  box  of  side  r  on  the  free 
storage  list.  If  the  free  storage  list  has  no  boxes  of  side  r,  then  let 
<\*  be  the  smallest  power  of  2  greater  than  r  for  which  the  free  storage 
list  has  a  box  of  side  q*.  (We  shall  show  that  when  ALLOCATE  is  called 
during  the  simulation,  q*  must  exist.)  For  q  ■  q*,q*/2, . . . ,4r,2r  in  order, 
select  a  position  x^  of  a  box  of  side  q,  delete  x^  from  the  list,  and 
add  to  the  list  positions  of  the  2d  pairwise  disjoint  boxes  of  side  q/2 
whose  union  is  C^.  Finally,  let  x  be  the  position  of  a  box  P  of  side  r 
on  the  free  storage  list,  delete  x  from  the  list,  and  return  the  value  x. 

The  time  taken  by  ALLOCATE  is  O(log^N) . 

Lemma  3.  Let  r  be  a  power  of  2.  Suppose  the  total  volume  of  boxes 
on  the  free  storage  list  is  at  least  r*  in  a  configuration  of  D  at  the 
beginning  of  a  call  to  ALLOCATE.  Then  this  call  can  produce  the  position 
of  a  blank  box  of  side  r  in  the  mass  store. 


Using  Leona  3,  which  can  be  proved  routinely  [L],  we  verify  chat 
ALLOCATE  can  always  find  the  requested  boxes  in  the  mass  store.  At  any 
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point  during  the  sinulation,  let  P^^,...  be  the  active  pages  in  the  mass 
store  and  m^.mg,...  be  the  values  of  their  nonblank  cell  counters.  Let  Pj 
be  assigned  to  block  5^.  According  to  Leona  2(a),  the  side  of  is  at 
most  tt ( (njj  log  N)^d).  The  mass  store  holds  smaller  inactive  pages  that 
were  assigned  to  in  previous  configurations  of  D.  Because  the  volumes 
of  these  smaller  pages  are  distinct  powers  of  2,  the  total  volume  of  pages 
ever  assigned  to  B^  is  at  most  twice  the  volume  of  P^,  hence  at  most 
2(TT((m^  log  N)^d))d.  By  Lemma  2(c),  die  tocal  volume  of  all  pages  in  the 
mass  store  is  at  most 


£ j2(TT((nij  log  N)1/d))dS  2(2d3*h  N  log  N)  <  (u*)d  . 

Ergo,  whenever  ALLOCATE  is  invoked  to  find  a  blank  box  of  side  r,  the  total 
volume  of  blank  boxes  on  the  free  storage  list  is  at  least  rd. 
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4.  RESEARCH  DIRECTIONS 

1.  Kosaraju  [Ko]  proved  that  multihead  one-dimensional  machines  with 
head-to-head  jumps  can  be  simulated  in  real  time  by  one -dimensional  machines 
with  just  one  access  head  per  worktape.  Do  multihead  tree  machines  and 
multidimensional  machines  with  head-to-head  jumps  share  this  property? 

Leong  and  Seiferas  [LS]  presented  real  time  simulations  of  nultihead 
multidimensional  machines  without  head-to-head  jumps. 

2.  How  fast  can  tree  machines  and  multidimensional  machines  simulate 
each  other?  Reischuk  [R]  devised  on-line  simulations  of  multidimensional 
machines  of  time  complexity  t(n)  by  tree  machines  in  time  0(t(n)c^°®*c^n^) 
for  a  constant  c.  Pippenger  and  Fischer  [PF]  gave  an  on-line  simulation 
of  a  tree  machine  by  a  one-dimensional  machine  in  time  0(t(n)  /log  t(n)). 

For  d  ^  2,  a  theorem  of  Grigoriev  [Gl  implies  that  every  tree  machine  can 
be  simulated  on-line  by  a  d-dimensional  machine  in  time  0(t(n)1+*^d 

but  no  better  upper  bound  is  known. 


J 
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